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Abstract. We have imaged the H92q (8.3 GHz), H75a (15 GHz), and H166« (1.4 GHz) Radio Recombination 
Lines (RRLs) from NGC 253 at resolutions of 4.5 pc (0.4") , 2.5 pc (0.2") and 53 pc (4.5") respectively. The 
H92q line arises from individual compact sources, most of which possess radio continuum counterparts. The line 
widths range from ~200 km s _1 for the sources near the radio nucleus to 70-100 km s _1 for the extranuclear 
ones. These lines are emitted by gas at a density ~10 4 cm -3 . The remainder of the cm-wave RRLs arise in lower 
density gas (~500 cm -3 ) with a higher area filling factor and with ten times higher mass. A third component of 
higher density gas (>10 4 cm -3 ) is required to explain the mm-wave RRLs. 

Key words. Galaxies: individual: NGC 253 - Galaxies: ISM - Galaxies: starburst - Radio lines: galaxies - Radio 
lines: ISM 



1. Introduction 

NGC 253 is a nearby (D=2.5 Mpc; 1"=12 pc) highly 
inclined spiral gal axy, undergoing: a vigo rous starburst 
(L F i R =3xl0 10 L^: lTelesco fc Harpedll980|l in the central 
^100 pc region. The ongoing star formation is concen- 
trated along an inclined ring. A number of compact near 
and mid-IR (NIR and MIR) sources have been discovered 
in this region (Forbes, Ward, and Depoy 1991; Forbes et al. 
1993; Pina et al. 1992; Keto et al. 1993; Sams et al. 1994). 
Optical imaging by Watson et al. (1996) revealed the pres- 
ence of four star clusters, which were identified with indi- 
vidual MIR/NIR sources (see also Forbes et al. 2000). The 
morphology of the radio continuum emission, which is un- 
affected by dust extinction, is quite different. Turner and 
Ho (1985) discovered a linear chain of compact sources at 
15 GHz, which were studied in detail by Antonucci and 
Ulvestad (1988) and Ulvestad and Antonucci (1997; here- 
after UA97). The latter were able to classify the brighter 
compact sources as synchrotron or free-free dominated, 
based on the derived spectral indices. Kalas and Wynn- 
Williams (1994) and Sams et al (1994) showed that most of 
the infrared sources are regions of low dust extinction and 
are not coincident with the radio sources. Hence, the ra- 
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dio sources and the infrared-optical sources seem to trace 
very different types of objects. In particular, the peak of 
the NIR emission is offset by about 3.5" southwest from 
the peak of the radio continuum emission (the radio nu- 
cleus). There is no associated radio emission towards the 
NIR peak and only weak NIR emission is observed near 
the radio nucleus. High resolution radio continuum ob- 
servations with the Very Large Array (VLA) suggest that 
the radio nucleus h osts an Active Galactic Nucleus (AGN; 
iTurner fc Holll985l UA97) , which was subsequ ently con- 
firmed l|Mohan et alJEool IWeaver et al.ll2002j) . The NIR 
continuum peak, identified as a region of active star forma- 
tion, also coincides with the position of the emission peaks 
of various ionised gas tracers ([Nell]: Baker, Krabbe and 
Storey 1998; Keto et al. 1999; Br7 : Forbes et al. 1993; 
optical continuum and line emission: Engelbracht et al. 
1998; Watson et al. 1996; Forbes et al. 2000). 

The immediate product of the starburst process, with 
its associated high star formation rate and efficiency, is 
a young massive stellar population. This population can 
studied through its surrounding ionised gas, the proper- 
ties of which are determined partly by the ionising stars 
themselves, and partly by gas dynamics. Recent star for- 
mation is well traced by NIR emission, except in regions 
of high dust extinction. However, such highly obscured 
star formation can be probed by radio continuum and 
hydrogen radio recombination line (RRL) emission (ex- 
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cept in cases where the continuum optical depth due 
to free-free emission is much larger than unity). RRL 
emission is insensitive to dust extinction, arises only in 
thermal gas, and can be imaged at arcsecond resolu- 
tion. NGC 253 was the first extragalactic object (ex- 
cludin g the Magellanic Cloud s) to have been detected in 
RRLs Ipeaauist fc Belli Il977fl . This line emission was ini- 
tially modelled a s being stimulated by the background 
nuclear emission Eh aver. Churchwell fc Walmslevt 11978: 
iMebold et all Il980^ . The first interferometric observa- 
tions of RRLs from NGC 2 53 were carried out by 
lAnantharamaiah fc Gossl l|l990l ) at higher frequencies (the 
H166a, HllOa and the H92a lines at 1.4 GHz, 4.8 GHz 
and 8.3 GHz respectively) using the VLA. Models indi- 
cated that these lines arise due to internal emission from 
high den sity gas (n P ~ 1 4 cm " 3 ; similar to compact HII 
regions). IPuxlev et alJ l|l997h observed mm- wave RRLs 
from this galaxy and suggested that the origin of the RRL 
emission at all frequencies is spontaneous emission from 
optically thin HII regions. A total ionisation rate of (3.7 ± 
0.8) xlO 53 s _1 was derived, an order of magnitude higher 
than that derived from infrared emission lines. The ioni- 
sation budget and hence the star formation rate as deter- 
mined by RRLs depends on the nature of the line emis- 
sion, which is not yet fully understood. High resolution 
observations (~1.3" or 16 pc) of the 8.3 GHz RRL from 
NGC 253 were carried out by lAnantharamaiah fe Gossl 
l|l996h . They showed that the line emission peaks at the 
radio nucleus (as do the higher frequency RRLs as well; 
IZhao et al.ll200l|) . They detected much weaker line emis- 
sion near the NIR peak, which is unexpected, since the 
infrared peak is probably the most active star forming 
site in this galaxy. Using the high resolution H92a and 
H75a data presented in this work, the RRL emission from 
the radio nucleus alone was studied separat ely and the re- 
sults were presented bv lMohan et alJl2002i Based on the 
results of the line emission models, it was concluded that 
the ionising source for the RRL emitting gas is probably 
a central low luminosity AGN with weak X-ray emission, 
which was su bsequently confirme d through CHANDRA 
observations llWeaver et al.U2f)02t . 

In order to investigate the nature of RRL emission in 
this galaxy, and hence the nature of obscured star forma- 
tion, we have performed high resolution observations of 
RRLs in NGC 253 at multiple frequencies using the VLA. 
The 8.3 GHz H92a line and the 15 GHz H75a line were 
observed using the VLA in the A configuration, which con- 
stitute the highest resolution RRL observations to date 
(0.4" and 0.2" respectively). These line emissions from 
the radio nucleus w ere studied separate ly and have been 
described in detail in lMohan et alJ 12002). The current pa- 
per deals with the remainder of this cm-wave RRL emis- 
sion. We have carried out these observations in order to 
(1) locate the sites of high frequency RRL emission and 
identify radio continuum counterparts to the line emitting 
sources, (2) to investigate the nature of the radio nucleus 
and the IR peak, based on their RRL emission, and (3) 



separate the diffuse and compact RRL emission in an at- 
tempt to construct a multi-density model of the ionised 
gas in NGC 253, using multi-frequency RRL data, sim- 
ilar to Arp 220 (jAnantharamaiah etall l2000(l . The 1.4 
GHz H166a line emission was observed at a resolution of 
~4" using the VLA in the Bn A configuration and the C 
config uration data presented in lAnantharamaiah fc Gossl 
(199QJ) were re-analysed. The aim was to better localise 
the large-scale line emission and hence (1) constrain the 
projected lateral size of the line emitting gas and thus 
obtain more accurate models, (2) investigate the relative 
importance of externally stimulated emission versus inter- 
nal (spontaneous) emission at 1.4 GHz compared to higher 
frequencies and (3) check the reality of the observed posi- 
tion offset between the l ow frequency and high frequency 
RRL emission found by lAnantharamaiah fc Gossl l|l990l) 
using the C and D configuration data. 

2. Observations 

2.1. Centimetre-wave RRLs: 8.3 GHz and 15 GHz 
observations 

The 8.3 GHz H92a and the 15 GHz H75a recombination 
lines from NGC 253 were observed using the VLA in the A 
configuration and the observational parameters are listed 
in Table ^ Only part of the 8.3 GHz dataset was used 
to study the weak line emission since the remaining data 
were affected by t he 3 MHz ripple in the VLA waveguide 
system 1 fhowever. lMohan et al.ll2002l used all of the data 
since the much stronger nuclear line emission is far less 
affected by the ripple). All data analysis were made using 
standard algorithms available within the software AIPS. 
The continuum emission at 8.3 GHz and 15 GHz was im- 
aged. Offline Hanning smoothing was appli ed to the line 
data to reduce the effects of Gibb's ringing l|Tavlor et all 
1999J). The line emission at both frequencies were imaged 
after applying corrections for the observed positional off- 
sets (see Sect. 13.1(1 . 

2.2. Decimetre-wave RRLs: 1.4 GHz data 

NGC 253 was observed at 1.4 GHz using the VLA in the 
BnA configuration and the observational details are listed 
in Table ^ The archival C configuration data were re- 
processed and are described in the appendix. The BnA 
array observations were performed in the 4 IF mode of 
the correlator in order to observe both the H166a and 
the H168a lines simultaneously. The continuum emission 
was imaged and deconvolved simultaneously with emis- 
sion from strong sources within the primary beam of 30'. 
The continuum emission from NGC 253 is surrounded by 
negative bowls of emission due to missing short spacing 
visibilities. Bandpass calibration every half hour proved 
adequate to track the time variation of the bandshape to 
an accuracy of <1 %. Both the H166a and the H168a 

1 See VLA Test Memorandum 158 by C.L. Carilli. 
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Parameter 


1.4 GHz 


8.3 GHz 


15 GHz 


VLA configuration 


BnA 


A 


A 


Date of observation 


01 Oct 95 


9 & 12 Jul, 


26 Jun 99 






9 Oct 99 




v TeBt of RRL (MHz) 


1424.7(H166a) 


8309.4 


15281.5 




1374.6(H168«) 


(H92a) 


(H75q) 


Bandwidth (MHz) 


3.03 


24.2 


46.9 


# of channels, # IF 


31, 4 


31, 1 


15, 1 


Spectral resolution" (km s _1 ) 


41.1 


56.4 


122.6 


Shortest baseline (kA) 


2.4 


7 


12 


Synthesized beam ("x") 


5.07x3.79 


0.50x0.28 


0.31x0.14 


Phase calibrator 


0023-263 


0116-219 


0118-272 


Bandpass calibrator 


0023-263,3048 


2251+158 


2251+158 


Peak continuum (mjy/beam) 


675 


40 


37 


Continuum noise (mjy/beam) 


0.13 


0.04 


0.13 


Noise per channel (mjy/beam) 


0.26 


0.12 


0.35 



The spectral resolution after off-line Hanning smoothing. 



lines were detected in the line cubes and the line param- 
eters were found to be consistent in the two datasets. 

3. Results 

3.1. Centimetre-wave RRLs: 8.3 GHz and 15 GHz 
observations 
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Fig. 1. 8.3 GHz continuum emission from the central disk 
of NGC 253 using VLA A configuration data. The rms 
in the image is 0.04 mjy/beam and the peak flux den- 
sity is 63 mjy/beam. The contour levels are (-5, 5, and 
higher in steps of 1.4) times the rms. The synthesized 
beam is 0.50" x 0.28", P.A.=0° (6 pc x3.4 pc). The sources 
5.49—42.3, 5.59—41.6 and the 15 pc region marked in the 
figure refer to H92a line emitting regions, discussed in 
Sect. The two remaining line sources 5.72—40.1 and 
5.75—39.9 cannot be seen clearly in the figure. 



The continuum images made at both 8.3 GHz and 15 
GHz are consistent with the images published by UA97; 
the continuum flux densities of individual compact sources 
agree to within 10 %. However, the positions of these 
sources in the 8.3 GHz and 15 GHz continuum images do 
not coincide. This positional discrepancy can be modelled 
as a relative global offset, indicating the uncertainty in 
the positions of the phase calibrators as a probable cause. 
Adopting the positions of compact sources in the A config- 
uration 15 GHz image of UA97, the positional offsets were 
derived to be +0.036"±0.01" in R.A. and +0.14"±0.01" 
in Dec for the 8.3 GHz data, and +0.021"±0.007" in R.A. 
and +0.012"±0.006" in Dec for the 15 GHz data, which 
were then applied to the data. All positions quoted in this 
paper are relative to 5.75— 41. 8 2 . The images of continuum 
emission at the two frequencies are shown in Fig. ^ an d 
Fig.0 

The integrated H92cv emission image using the full res- 
olution dataset is shown in Fig. |3Ji overlaid on an image 
of the continuum emission from the region. The emission 
from the inner 25 pc region is shown in Fig. |3p, over- 
laid on a higher resolution 15 GHz continuum image. The 
line emission is spatially extended, with the peak of the 
emission coinciding with the peak of the radio contin- 
uum (radio nucle us). This nuclear lin e emission was dis- 
cussed in detail in lMohan et alJ (|2002^ . The strongest line 
emission arises in an extended region surrounding the nu- 
cleus, spread over 15 pc. In addition, multiple compact 
line emitting sources were found to lie along the larger 
continuum disk. The morphology of the H92ev line emis- 
sion is quite complex and is discussed in detail in Sect.^ 
The velocity field of the ionised gas in the 15 pc region 
indicates that the velocity gradient is almost in the east- 
west direction and the value increased from 200 to 230 

2 The corrected (B1950.0) coordinates of 5.75-41.8 are 
00 h 45 m 05 8 .751, -25°33'41".85 at 8.3 GHz and 00 h 45 m 05 s .750; 
-25°33'41".86 at 15 GHz. 
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Fig. 2. 15 GHz continuum emission from the central disk 
of NGC 253 using VLA A configuration data. The rms 
in the image is 0.12 mJy/beam and the peak flux density 
is 44 mJy/beam. The contour levels are (-3, 3, 5, and 
higher in steps of 1.35) times the rms. The synthesized 
beam is 0.31"x0.14", P.A.=-20° (3.7 pc xl.8 pc). The 
sources 5.49—42.3, 5.59—41.6 and the 15 pc region marked 
in the figure refer to H92a line emitting regions, discussed 
in Sect. 0J The two remaining line sources 5.72—40.1 and 
5.75—39.9 cannot be seen clearly in the figure. 



3.2. Decimetre-wave RRLs: 1.4 GHz data 

The 1.4 GHz continuum image is shown in Fig. [SJ The 
data for both the RRLs were combined and the line emis- 
sion was imaged after offline Hanning smoothing. The line 
emission coincides with the position of the peak of the 
continuum emission and is marginally spatially resolved. 
The line does not increase in strength by more than 10 % 
in images made at lower resolution. Hence the synthesized 
beam (5.1"x3.8") covers almost all of the line emitting re- 
gion. The spectrum obtained from the combined dataset 
towards the continuum peak is shown in Fig. EJ The line 
parameters are listed in Table 
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km s" 1 from east to west, consi s tent w ith the observa- 3 
tions of I Anantharamaiah fc Goss] {l996). The analysis of ^ 
the velocity field is not undertaken in this paper. 

The 15 GHz line cube was used to search for sources of 
H75a line emission. The channel images of the H75a line 
emission as well the H75a spectra towards the two sources 
are shown in Fig. 0] The peak of the line emission is coin- 
cident with the radio nucleus and the only other source of 
RRL emission arises from a compact source ~0.35" north- 
west of the nucleus. This source seems to have correspond- 
ing H92a line emission at 8.3 GHz, associated with the 
continuum source 5.76—38.9 (Sect. 0J), and is the second 
strongest line emitting source at the latter frequency. 
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Fig. 3. The four extranuclear H92a sources, discussed 
in Sect. 14.11 and the central extended emission are de- 
tected. Top: The H92a line emission (thick contours) 
overlayed on the 8.3 GHz continuum emission (thin con- 
tours). Two compact sources of line emission (5.49—42.3 
and 5.59—41.6) are seen towards the SW of the emis- 
sion surrounding the radio nucleus. The synthesized beam 
is 0.50"x0.28", P.A.=0° and rms of the continuum im- 
age is 0.04 mJy/beam. Bottom: A close-up view of the 
H92a line emission towards the 25 pc region (grey scale) 
overlayed on a high resolution 8.3 GHz continuum im- 
age (contours). The line emission consists of two addi- 
tional extranuclear sources (5.72—40.1 and 5.75—39.9) and 
a contiguous 15 pc sized central line emitting region. 
The synthesized beam of the continuum and line images 
are 0.35"x0.19", P.A.=-3° and 0.45"x0.24", P.A.=-3° 
respectively. The rms of the continuum image is 0.07 
mJy/beam and the contours are in steps of 1.4 starting 
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Source 


Peak flux 


Central velocity" FWHM" 


Integrated flux 




J • J f T /"I \ 

density (mjy/beamj 


(km s" 1 ) 


(km s" 1 ) 


(X10 W m ) 




i-4 Kjiiz rLiOocy. Line 


emission (BnA 


configuration) 




7.4" x 5.2" area 


1.8 ± 0.1 


200 ± 20 


160 ± 14 


1.6 ± 0.25 




8.3 GHz H92a line 


emission^ (A configuration) 




5.49-42.3 


0.5 ± 0.1 


201 ± 7 


73 ± 17 


1.2 ± 0.3 


5.59-41.6 


0.5 ± 0.1 


172 ± 7 


65 ± 16 


1.0 ± 0.3 


5.72-40.1 


0.46 ± 0.09 


170 ± 15 


101 ± 23 


1.4 ± 0.3 


5.75-39.9 


0.54 ± 0.09 


175 ± 10 


112 ± 20 


1.8 ± 0.3 


'15 pc region' 


3.2 ± 0.15 


194 ± 5 


211 ± 12 


20 ± 2 




15 GHz HI 5a line 


emission (A configuration) 




Nucleus 


1.8 ± 0.3 


224 ± 17 


230 ± 40 


23 ± 3 


5.76-38.9 


1.6 ± 0.3 


205 ± 18 


230 ± 43 


20 ± 3 



a Heliocentric velocity, optical definition. 

b Corrections for finite velocity resolution not applied. 

c All sources identified by coordinates are spatially unresolved. See Table^for synthesized 
beam sizes. 
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Fig. 4. (a). VLA A configuration 'naturally weighted' im- 
ages of H75a emission in individual channels from the 
central 15 pc region of NGC 253. The central heliocen- 
tric velocity (in km s _1 ) of each channel is marked in each 
frame. The synthesized beam is 0.35"x0.22", P.A.=-10° 
and the rms in the channel images is 0.36 mJy/beam. The 
contour levels are (3.5, 3.8, 4.1, 4.4, 4.7, 5, 5.3) times the 
rms. (b) . H75ev spectrum towards the continuum peak (ra- 
dio nucleus), (c). H75cv spectrum towards the line source 
5.76-38.9. 
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Fig. 5. 1.4 GHz continuum emission from the inner disk of 
NGC 253 using the VLA BnA configuration data. The rms 
in the image is 0.13 mJy/beam and the contour levels are 
from 1.3 to 664 mJy/beam in steps of 1.7. The synthesized 
beam is 5.1"x3.8", P.A.=51°. 



The RRL emission f r om N GC 253 was imaged by 
lAnantharamaiah fc Gossl l|l990|) at low resolution, using 
the VLA in the C and D configurations. The 1.4 GHz 
line emission was found to be shifted SE from the peak of 
the continuum and the 4.8 GHz and 8.3 GHz line emis- 
sions. Subsequently, the D configuration data were found 
to exhibit problems in the spectral domain and was hence 
unusable. The reanalysed C configuration data (appendix) 
shows the line emission to be centred on the continuum 
peak. Hence the position al offset for the 1.4 GHz lin e emis- 
sion presented earlier bv lAnantharamaiah fc Goss] {l990) 
is not confirmed. 
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HELIOCENTRIC VELOCITY (km/s) 

Fig. 6. Co-added H166a and H168a spectra towards the 
continuum peak of NGC 253 using the VLA BnA config- 
uration data for the natural weighted image. Synthesized 
beam is 5.1"x3.8", P.A.=50° and the rms noise in the line 
cube is 0.26 mJy/beam. 



4. The morphology and analysis of the extended 
H92a RRL emission 

In this section, the morphology of the extended H92a line 
emission is analysed in order to identify the individual line 
emitting sources and derive their properties. Continuum 
counterparts can be associated with most, but not all, of 
the line emitting sources. We identify each of the RRL 
sources by the continuum counterpart or the nearest con- 
tinuum source. There are two compact line sources to the 
SW (Fig. Eli), and the central region consists of two ad- 
ditional compact sources along with the extended 15 pc 
sized nuclear emission (Fig. 0t>) ■ These four extranuclear 
line sources are described in the next section, followed by 
a discussion of the emission from the central 15 pc region. 



(c) 5.72-40.1 
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Fig. 7. H92a VLA spectrum towards the four compact 
extranuclear line sources using the A configuration data, 
identified by the coordinates of the associated/nearest 
continuum source : (a). 5.49-42.3, (b). 5.59-41.6, (c). 
5.72-40.1 and (d). 5.75-39.9. The synthesized beam is 
0.50" x 0.28", P.A.=-2°. The rms in the line cubes is 0.12 
mJy/beam and the velocity resolution is 56.4 km s _1 . 

ated with the IR peak (see Sect.s l5T2l and l6.lfl . Two addi- 
tional sources of compact H92a line emission are detected 
in the central region (Fig.Qb). The western source is close 
to the continuum source 5.72—40.1. The latter exhibits a 
flat spectral index from 4.8 GHz up to 23 GHz and was 
identified by UA97 as the brightest thermal source in NGC 
253. The eastern source, however, has no obvious contin- 
uum counterpart, and in fact lies in a region of minimum 
continuum emission. Nevertheless, the nearest continuum 
source 5.75—39.9 will be used to refer to this line emitting 
source. The H92a spectra towards these extranuclear line 
sources are presented in Fig. [3 These line sources are spa- 
tially unresolved and the spectral parameters are listed in 
Table 



4.1. The extranuclear H92a sources 

The western most line source (seen in Fig. [SJi) is nearly 
coincident with the continuum source 5.49—42.3. A sig- 
nificant fraction of its continuum emission is thermal 
(the 8.3 GHz-23 GHz spectral index is -0.28 ± 0.36; 
UA97), supporting its association with the line emis- 
sion. The line emitting region 1.5" east of this source 
(Fig. Ot.) is displaced ~0.1" southwest of the compact con- 
tinuum source 5.59—41.6. This object is the second most 
prominent RRL source i n the l ower resolution image of 
lAnantharamaiah fc Gossl l|l996|) and is probably associ- 



4.2. The H92a emission in the central 15 pc region 

lAnantharamaiah k. Gossl l|l996f) found that the strongest 
H92a emission arises from the region surrounding the nu- 
cleus. Our observations, which have three times improved 
angular resolution, can spatially resolve this central emis- 
sion, which peaks at the radio nucleus. High resolution 
channel images of this 15 pc region are shown in Fig. [SJ 
The line emission in this region can be resolved into four 
discrete components, as can be seen in Fig. El though the 
identification of the individual sources is uncertain due 
to the complex morphology and the low signal to noise. 
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However, each of the sources can be best identified in dif- 
ferent channel images. One or more channels were identi- 
fied where the sources can be distinguished from the sur- 
rounding emission (the southern source at a velocity of 
119.0 and 90.8 km s _1 , the westernmost at 203.6 km s _1 
and the two northern sources at 231.9 and 260.1 km s" 1 ). 
These channels were used to derive the peak positions by 
fitting spatially unresolved gaussians using the AIPS task 
JMFIT. The errors in the positions were determined by 
the fitting procedure. Given the low signal-to-noise of the 
emission (the peak emission is in the range of 5-9 a), a si- 
multaneous 3-d deconvolution of all four line sources was 
not attempted. 

The derived positions of these four sources are listed 
in Table [3] and are also marked in Fig. which shows 
a high resolution image of the integrated line emission 
from this region. The offsets between the marked posi- 
tions and those of the line peaks in the figure are due to 
blending of emission from the four sources. In Fig. Et: the 
integrated line emission is overlaid on a high resolution 
15 GHz continuum image. Apart from the nucleus, the 
source immediately to the west coincides with a strong 
continuum source that is not listed in the catalogues of 
UA97. This H92a line source also has a H75a line emission 
counterpart (Fig.0J. The remaining two H92a sources do 
not have obvious continuum counterparts. The H92a spec- 
tra towards each of the four sources are shown in Fig. 1101 
and their line parameters are summarised in Table 0] 
The sum of the H92a emissions from the four sources is 
(18±1) x 10~ 23 W m~ 2 , comparable to the total emission 
of (20 ± 2) x 10~ 23 W m -2 from the entire region. Hence 
all of the extended emission observed within the 15 pc 
region can indeed be accounted for in terms of four indi- 
vidual compact sources. 

lAnantharamaiah fc Gossl l)l990|) imaged the H92a 
emission at a resolution of 7.4" x 5.2" and detected a total 
line emission of 44xl0~ 23 W m~ 2 , whereas our observa- 
tions account for only 25.5xl0 -23 W m -2 over the entire 
disk. The remainder must be spatially extended so as to 
fall below the surface brightness sensitivity of the current 
observations. In this section, we have decomposed the ob- 
served H92a emission into eight distinct compact sources. 
The positional offsets between these line sources and the 
associated continuum sources are typically <0.1". These 
offsets are therefore significant when compared to the size 
of the synthesized beam. The continuum emission in NGC 
253 is dominated by non-thermal radiation, whereas RRLs 
arise in thermal gas. Hence a spatial offset between the two 
emission regions could w ell exist, as in the case of NGC 
1808 at high resolution l(Kotilainen et all Il996h . These 
offsets are probably due to lower surface brightness of 
thermal compared to nonthermal emission, coupled with 
source confusion due to the high surface density of com- 
pact sources. Continuum emission imaged at a resolution 
higher than that of the current images should be able to 
delineate the actual continuum counterparts of the RRL 
sources. 
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00 45 05.84 05.80 05.76 05.72 
RIGHT ASCENSION (B1950) 



Fig. 9. (a). High resolution integrated VLA H92a line 
emission from the 15 pc region, summed over veloc- 
ities 316.6 to 90.8 km s _1 . The synthesized beam is 
0.38"x0.20", P.A.=-2°. The contour levels are (1, 1.4, 
1.8, 2.2, 2.6, 3. 3.4, 3.8, 4.2, 4.6, 5) times 10 3 Jy/beamxHz. 
The 'star' symbols mark the derived positions of the four 
line sources, as listed in Table (b). The positions of 
the four H92a sources are marked as 'star symbols' on 
the contours of 15 GHz continuum emission imaged at 
high resolution. The 15 GHz continuum image is aligned 
using the continuum and line emission from the radio nu- 
cleus. The synthesized beam of the continuum image is 
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Fig. 8. VLA A configuration images of the H92a line emission in individual spectral channels in the central '15 pc 
region' of NGC 253, at higher resolution. The central heliocentric velocity of each channel is marked in each frame. Line 
emission in every alternate channel is independent due to Hanning smoothing. The synthesized beam is 0.41" x 0.22", 
P.A.=10° and the rms in the channel images is 0.13 mJy/beam. The contour levels are (-3, 3, 3.5, 4, 4.5, 4.75, 5, 5.5, 
5.8, 6, 6.5, 7, 7.5) times the rms. 

Table 3. Positions of the H92a and H75a line emitting sources in the 15 pc region 



Source 


RA 


Dec 




(B1950.0) 


(B1950.0) 




The four H92a sources 


(see Fig. 0) 


1 


S .002 (2) 


2".90 (5) 


2 


s . 022 (2) 


2" .31 (5) 


3 


a .026 (1) 


2".87 (2) 


4 


s . 044 (1) 


2".80 (2) 




The second H75a 


source 


5.76-38.9 S .019 (2) 


-2".94 (5) 



NOTE: The coordinates are with respect to 00 h 45 m 05 H .751; 
-25°33'41".85, which is the (B1950.0) position of 5.75-41.8 
in the 8.3 GHz image after applying positional corrections 
fScct. l3~Tl . Positive offsets are for eastward and northward di- 
rections. The gaussian fitting errors in the final decimal place 
are given in parenthesis. 



5. A multi-density model for the ionised gas in 
NGC 253 

RRL data can be used, along with the associated radio 
continuum measurements, to model the physical proper- 
ties of the line emitting (photo-)ionised gas. The atomic 
level populations are derived 3 assuming they are not in 

3 The departure c oefficients are tak en from 
a program by ISalem fc Brocklehursd lll979l) . 



local thermodynamic equilibrium (LTE). The ionized gas 
is modelled either as a uniform sphere, a slab or a col- 
lection of spherical HII regions, all of constant density. 
Density, electron temperature, size (radius for spherical 
models and lateral size and thickness for slab models) and 
the number of HII regions (for the collection of HII re- 



modified bv IWalmslev fc Watsonl dl982l) and 
IPavne. Anantharamaian^^5rteksonT(^93) . 
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Table 4. Parameters for the H92a sources in the 15 pc region 



Source 
name 




Peak line 
flux density 
(mjy/beam) 


Central 
velocity" 
(km s _1 ) 


FWHM" 

(km s" 1 ) 


Integrated flux 
(xl0~ 23 W m" 2 ) 
H92a H75a 


Continuum 2 
(mjy/beam) 
8.3 GHz 15 GHz 


Source d 
size 


source 


r 


0.5 ± 0.1 


-220 


~150 


~2 


<9.1 


10.5 


9.5 


unres. 


source 


2 


0.6 ± 0.1 


180 ± 10 


220 ± 20 


3.9 ± 0.4 


<9.1 


30 


26 


unres. 


source 


3 


0.75 ± 0.06 


200 ± 10 


240 ± 22 


5.1 ± 0.9 


20 ± 3 


34 


35 


unres. 


source 


4 


1.15 ± 0.06 


192 ± 6 


198 ± 13 


6.7 ± 0.5 


23 ± 3 


64 


59 


unres. 



Heliocentric velocity, optical definition. 

Corrections for finite velocity resolution not applied. 

Continuum flux per synthesized beam measured at the position of the line emission. 

Spatially unresolved, angular size <0.35". 

The H92a parameters of source 1 are less certain. 




600 400 200 -200 600 400 200 -200 
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-0.2 u ' 1 ' 1 ' 1 ' 1 1 -0.2 L- 1 - 
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Heliocentric Velocity (km/s) Heliocentric Velocity (km/s) 

Fig. 10. The VLA A configuration H92a spectra towards 
the line emitting sources in the 15 pc region, obtained 
towards the positions listed in Table |21 The synthesized 
beam is 0.51" x 0.28", P.A.=-2°. The rms noise in the 
channel images is 0.12 mjy/beam. 



gions model) are free parameters. The temperature of the 
photo-ionized gas in all the models is taken to be 7500 
K, a value typical of galactic HII regions, since the re- 
sults of the models do not depend significantly on the 
temperature. The values for the continuum and line flux 
densities are calculated for a grid of values of the free 
parameters and the observed flux densities are used to 
constrain the solution space. The range of derived val- 
ues of the free parameters is determined primarily by 
the observational errors in the flux densities. Externally 



stimulated emission is also consid ered (in addition to in- 
terna lly stimulated emission, see lAnantharamaiah et alJ 
1993) using the observed value of the continuum flux den- 
sity. It is found that the derived range of the free param- 
eters allow for substantial externally stimulated emission 
as well, but the values of these parameters are not af- 
fected significantly due to limits to the continuum bright- 
ness temperature. Further det ails of the modelling pro- 
cedure have been described in iMohan et alJ (2001), and 
references therein. In Sect. 15.11 we attempt to model the 
entire low-resolution multi-frequency RRL data as aris- 
ing from a single component of ionised gas with uniform 
properties. In Sect. 15.21 we use our higher resolution ob- 
servations of cm-wave RRLs to derive the properties of 
the individual line emitting regions, and finally, in Sect. 
4.3, we propose a multi-density model for the ionised gas 
component in NGC 253. 

5.1. Modelling the low-resolution RRL data 

IPuxlev et all l|l997f) detected the 99 GHz H40a line in 
NGC 253 using single-dish observations and also com- 
piled low resolution interferom etric RRL observations at 
1.4 G Hz, 4.8 GHz and 8.3 GHz l|Anantharamaiah fc g"ossI 
Il990|) . They show that the RRL data from 1.4 GHz up to 
99 GHz is well fit by a F^ ~ v 2 curve (that is, the line 
spectral index a— 2) . In fact, our revised estimate of the 
1.4 GHz line flux density lies much closer to the best fit 
line than the earlier measurement. These data are plotted 
in Fig.^] A function of the form F u oc v a is fit to the data 
and the best fit value of a is 1.93 ± 0.04. The theoreti- 
cally expected dependence of F„ on frequency when t c <S^1, 
| T; | 1 , assuming LTE conditions and negle cting pressure 
broad ening, is indeed F„ ~ v 2 , as noted bv lPuxlev et alJ 
lll997h . They therefore suggest that all of the RRL emis- 
sion from 1.4 GHz up to 99 GHz arises in LTE gas which 
is optically thin. This result is surprising since we do not 
expect either LTE conditions or optically thin conditions 
to hold for ionised gas over this range of frequencies. 

In Fig. El we show the theoretical contours of the spec- 
tral index a between 1.4 GHz and 99 GHz as a function 
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Fig. 11. Integrated line flux density vs. frequency for 
NGC 253 for low resolution RRL data. The error bars 
are marked and the solid line is the best linear fit through 
the da ta. 1.4 GHz data : our work; 4.8 GHz and 8.3 GHz 
data : | A nantha ramaiah fc Gossl {l990); 99 GHz data : 
IPuxlev et alJ l)l997|) . 



of electron density and emission measure EM, including 
pressure broadening, for LTE as well as for non-LTE con- 
ditions. In the LTE case, a is quite close to 2.0 for EM<10 6 
pc cm~ 6 (or t c at 1.4 GHz <0.4) and increases for higher 
values of EM due to optical depth effects. The weak de- 
pendence of a on density for low values of EM is due to 
pressure broadening. However, the behaviour of a is very 
different for the non-LTE case. The value of a is closer to 
2.0 for certain high-EM regions in the density-EM plane 
compared to the LTE case. The ratio between the non- 
LTE and LTE flux densities is l/b n (b n is the departure 
coefficient for level n), and is hence close to unity, only 
for optically thin regimes. He nce, an LTE analysi s of the 
line emission, as performed bv lPuxlev et alJ l|l997f) . based 
on the observed v 2 dependence, will not be valid if the 
optical depth is >0.5. 

Therefore we investigate whether the observed line flux 
densities and the v 2 dependence of the integrated RRL 
flux from NGC 253 can indeed be modelled as emission 
arising from a single parcel of gas. The ionised gas is mod- 
elled as a rectangular slab at T e =5000 K, with the density 




-2 -1 1 2 3 4 5 
Density log(cm -3 ) 



1 2 3 4 5 6 
Density logfcm -3 ) 



Fig. 12. The expected dependence of the line spectral in- 
dex a (defined as Y v oc v a , where F„ is the integrated 
flux density of the RRL at frequency v) between the 1.4 
GHz H166a and the 99 GHz H40a line flux densities as a 
function of density and emission measure under (a) . LTE 
conditions and (b). non-LTE conditions. The boxes mark 
the solution space for Components I (open box) and II 
(hatched box) of the ionised gas, derived in Sects 15.21 and 



n e , lateral extent and thickness being the free parameters. 
The gas is required to produce the observed integrated 
line flux densities at the four frequencies. The valid solu- 
tions obtained violate two observational constraints. The 
minimum lateral angular extent (^10") is larger than the 
observed value (<6"). Additionally, the predicted ther- 
mal continuum at 99 GHz is much hi gher than t he es- 
tim ated value from th e observations of ICarlstroml l|l990|) 
and IPeng et ah (199G). In a model where the calculated 
flux densities are normalised, and only the calculated line 
spectral index a is constrained to be consistent with the 
observed value of a=1.93 ± 0.04, no solutions are ob- 
tained for LTE emission within 2<r. The acceptable non- 
LTE models in this case are confined to n e <100 cm" 3 and 
EM<10 4 with some isolated high density solutions as well. 
An observed spectral index of two, therefore, cannot be 
used to infer emission measures assuming LTE conditions. 

Hence, all of the RRL emission from 1.4 GHz up to 
99 GHz cannot arise in the same component of ionised 
gas. However, the morphology of the line emission in the 
H92a line images is observed to be complex and consists 
of individual sources, which may have different properties 
(for example, the line from the radio nucleus arises in gas 
at a density 10 times higher than derived in the above 
models and contributes abo ut 10 % of the total ionisation 
derived: iMohan et aLlEool^K Therefore each of the com- 
pact line emitting sources must be modelled separately, 
as is done in the following section. 
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5.2. Modelling the cm-wave RRLs: High density 
ionised gas 

The physical properties of the compact H92a emitting 
sources can be modelled in a more realistic fashion with 
our high resolution data since the smaller synthesized 
beam size strongly constrains the projected lateral size 
of the ionised gas. The four extranuclear line sources and 
the four sources inside the central 15 pc region are mod- 
elled individually. The models are constrained to produce 
the observed H92a and H75a flux densities. The free- 
free emission at 8.3 GHz and 15 GHz is also constrained 
to be less than the measured continuum flux densities. 
Further, the angular size of the gas in these spatially unre- 
solved sources is constrained to be <0.35", corresponding 
to the size of the synthesized beam. The observational con- 
straints used in the modelling are summarised in Table 
for the extranuclear sources and in Table 01 for the sources 
in the 15 pc region. The valid solutions obtained for the 
extranuclear sources are listed in Table [5] and notes on 
each object are given below. 

(1) 5.49—42.3: No solutions are obtained for a spherical 
HII region model; uniform slab models yield solutions cor- 
responding to thi n slabs. Since this s ource is not listed as 
a Ha emitter bv iForbes et al.1 l|2000|) . the integrated Ha 
flux of 2xl0 -15 W m -2 derived from the ionisation rate 
implies a (foreground) extinction of Ay > 10. 

(2) 5.59-41.6 ( = IR peak ?): \Vmn et alJ lll992h . 
ISams et alJ l|l994h and others have shown that the ra- 
dio nucleus and the IR peak are two separate objects, 
and that astrometric errors do not allow an unambiguous 
association of any other radio source with the IR peak. 
The closest radio sources are 5.59—41.6 and 5.62—41.3, 
which have spectral indices of +0.61 ± 0.15 and —0.21 
± 0.14 respectively, and are hence dominated by thermal 
gas. However, the detection of recombination lines from 
5.59—41.6 leads us to suggest that this is indeed the IR 
peak. This issue is discussed further in Sect. 16. ll 

(3) 5.72—40.1: We model the line emission from this ther- 
mal object along with the continuum constraints derived 
by UA97. The derived gas mass, density, and size are simi- 
lar to the UA97 values. However, only an insignificant frac- 
tion of the continuum at j/<8.3 GHz is thermal, whereas 
the continuum spectral index is observed to be flat at fre- 
quencies > 4.8 GHz. Hence, a single-density model cannot 
explain the observed flat spectral index as well as the RRL 
emission. Using the estimated Ha strength of 1.2 xlO -14 
W m~ 2 and ident ifying the source with a Ha knot in 
IForbes et al.l (|2DO0h ~0.4" eastwards, the (foreground) ex- 
tinction Ay is estimated to be 13. 

(4) 5. 75— 39. 9: The line emitting gas is modelled assuming 
that the free-free flux density of the thermal gas does not 
exceed the measured continuum of 5.5 mJy at 8.3 GHz. 
No Ha or IR emiss ion knots are observed in this region 
(jForhes et all l200nft and hence the extinction Ay must be 
> 9. 



The solutions obtained for the sources in the central 
15 pc region are similar for a uniform sphere as well as a 
rectangular slab geometry, and are described in Table El 
The derived range of densities for sources 1 and 2 are wider 
than those for sources 3 and 4 since only upper limits to 
the H75a line flux density are available for the former two. 
The spectral indices of RRL flux density between pairs 
of frequencies from 1.4 GHz to 99 GHz derived from the 
models range between 0.5 and 5.0 which differ significantly 
from 2.0 (see Sect 15 .If) . The derived properties of the four 
sources (Table [SJl are similar. 

The gas densities of the eight compact H92a sources lie 
in the range 5xl0 3 -10 4 cm -3 , similar to compact HII re- 
gions in our galaxy. To illustrate the advantage of high res- 
olution observations (which set stricter constraints on the 
projected size), we also model the 15 pc region as a whole. 
The results obtained are different compared to those de- 
rived above (lower densities, larger sizes and a dominant 
contribution from stimulated emissi on). However, t he re- 
sults are similar to those derived bv IPuxlev et al.1 l|l997t) 
from low resolution data. Hence, high resolution observa- 
tions are essential for reliably determining the properties 
of the compact component of ionised gas. 

5.3. Modelling the decimetre-wave RRLs: multi-density 
medium 

The sum of the observed continuum and line flux densities 
of the eight H92a sources, derived in Sect. 15.21 are listed 
in Table[7|(as Component I). This component contributes 
no more than 2 % (or 21 % for the low density models) 
of the observed H166a emission at 1.4 GHz, between 18- 
36 % of the HI 10a emission at 4.8 GHz and about 9 % 
of the H40a emission at 99 GHz. Therefore, in order to 
explain the remaining RRL emission at these frequencies, 
additional components of ionised gas, different from that 
detected in our high resolution H92a data, must exist. 
No plausible models can account for the remainder of the 
1.4 GHz as well as the 99 GHz line emissions to arise in 
the same component of gas. Ionised gas at densities >10 4 
cm" 3 is needed to account for the observed mm- RRL emis- 
sion, and this gas would produce negligible line emission 
at 1.4 GHz due to continuum optical depth effects. Hence, 
a minimum of three components of ionised gas seem to 
be needed to explain all of the RRL data on NGC 253 
: Component I which produces the high brightness part 
of the H92a emission, described in the previous section, 
Component II which explains most of the 1.4 GHz line 
and the remainder of the 4.8 GHz and 8.3 GHz lines, and 
Component III, which leads to the mm- wave line emission. 

Component II is modelled to produce all of the remain- 
ing 1.4 GHz, 4.8 GHz and 8.3 GHz line emissions, with the 
99 GHz line flux density imposed as an upper limit. The 
transverse size is constrained to be <7.5", the size of the 
low resolution emissions. The derived electron densities 
are between 300-1000 cm -3 , effective sizes between 30- 
60 pc and ionising photon rate is (3-6) xlO 52 s -1 . The 
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Table 5. Model results for the extranuclear H92cv sources 



Parameter 


5.49-42.3 


5.59-41.6 


5.72-40.1 


5.75-39.9 


Density (xW cm _J ) 


4-10 


2-20 


2-10 


2-10 


Effective size (pc) 


2-3. 5 a 


1-5 


1.4-2.7 


1.7-2.9 


r c at 8.3 GHz 


0.2-1.5 


0.2-3 


0.2-1.2 


0.1-1.5 


S ff at 8.3 GHz (mjy) 


0.5-0.7 


1-2.5 


2-3 


2-3 


N Lyc (xlO 51 ph. s" 1 ) 


1 


2 


1-10 


2-5 


Mass (xlO 3 M Q ) 


0.16-0.36 


0.3-1 


0.35-1.3 


2 


4.8 GHz free-free flux density (mjy) 


0.2-0.4 


0.7-2 


1-4 


1-2.5 


99 GHz free-free flux density (mjy) 


0.7-0.9 


1.3-2 


2-3.5 


2-4 


1.4 GHz line flux ratio 6 (%)' 
4.8 GHz line flux ratio 6 (%) 





0.1 


<1 


0.05 


0.2-0.8 


0.7-1.8 


1-4 


1-2.5 


99 GHz line flux ratio 6 (%) 


0.2-0.3 


0.4-0.7 


0.7 


0.7 



Only thin slab solutions with thickness between 0.01-0.04 pc. 

Ratio of predicted RRL flux of the eight H92a sources to the observed value. 



Table 6. Model results for the H92a sources in the 15 pc region 



Parameter 


Source l a 


Source T 


Source 3 


Source 4 




Low 


High 


Low 


High 






Density (xl0 a cm -3 ) 


2 


10 


0.9 


6 


8-10 


7-10 


Effective size (pc) 


5 


1.8 


6.4 


2.6-3.7 


2-2.5 


2.5 


t c at 8.3 GHz 


0.2 


1.7 


0.05 


0.8 


1-2.5 


1-2 


Sff at 8.3 GHz (mjy) 


0.4 


0.3 


2 


4 


6 


5.5 


N Lyc (xlO 51 ph. s" 1 ) 


2.5-5 


4-5 


1.5 


4 


8-15 


6-14 


Mass (xlO 3 M Q ) 


3-5 


0.8 


3.2 


1.3 


1.5-2.5 


1.3-2 


4.8 GHz free-free flux density (mjy) 


3-5 


1.6 


1.2 


3-5 


2.7 


2.5-5 


99 GHz free-free flux density (mjy) 


2.3-4.3 


4 


1.8 


2-7 


6-10 


6-12 


1.4 GHz line flux ratio 6 (%) 
4.8 GHz line flux ratio 6 (%) 


1 


0.05 


19 


0.1 


0.05 


0.05 


5 


1.6 


13 


5 


3.5 


5 


99 GHz line flux ratio 6 (%) 


0.8 


1.1 


0.3 


1 


2.5 


2.5 



a For sources 1 and 2, the density range is large; typical low and high density solutions are listed. 
6 Ratio of predicted RRL flux of the eight H92a sources to the observed value. 



stimulated emission component due to the background 
nonthermal radiation is between 50-90 %. If the stimulated 
emission is set to zero, we obtain solutions for densities 
between 200-500 cm~ 3 , but with very similar ionisation 
rates. The derived properties of the gas in Components I 
and II are listed in Table[7] This gas (Component II) con- 
tributes 10-12 % of the remaining 99 GHz emission. Hence 
a third higher density component is obviously necessary . 
Since the H40a line was detected bv lPuxlev et al 1 lll997t) 
at low signal to noise, we do not attempt to model this 
high frequency line emission in detail. 

6. Discussion 

Though high resolution radio continuum observations of 
NGC 253 exists, recombination lines provide additional 
advantages. The local gas density can be derived from 
RRLs whereas continuum data yield only the rms den- 
sity, and the densities derived in this work are higher 
than the rms densities for some of the thermal sources. 
Moreover, the thermal fraction of the continuum emission 



is not known accurately due to the large errors in the ob- 
served continuum spectral indices. Such errors also result 
to difficulties in determining free-free optical depth effects, 
especially in the case of multiple density components in- 
side a single source. 

The continuum spectral index measurements by UA97 
indicate that thermal gas is present over a substantial area 
and also forms a significant part of the continuum emis- 
sion from a number of compact sources. Hence, it is no- 
table that 8.3 GHz and 15 GHz line emission is detected 
from only a handful of sources, which may be ascribed 
partly to lack of adequate sensitivity. Additionally, since 
the observed RRLs are sensitive primarily to gas at den- 
sities 10 3 -10 4 cm -3 , the line emission would not be traced 
efficiently by high resolution observations of the H92a line 
emission if a substantial part of the thermal gas detected 
in the continuum is at very different densities. Lower den- 
sity diffuse gas, for example, can give rise to cm- wave RRL 
emission which is below the sensitivity limit of our obser- 
vations while still being able to make the observed contin- 
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Observational parameter 


Component I 


Observations 
(low res.) 


Comp.s II+III 
(expected) 


H166a at 1.4 GHz (xlO -23 W m~ 2 ) 


0.03-0.3 


1.6 


1.5 ± 0.4 


TT1 1 n j AO /"I TT / — 23 TT7 — 2 \ 

HllOa at 4.8 GHz (xlO W m ) 


3.4-6.9 


19.1 


14.5 ± 4 


H92q at 8.3 GHz (xlO~ 23 W m~ 2 ) 


25.5 


44.4 


19 ± 3 


H40q at 99 GHz (xl(T 23 W m" 2 ) 


564 


6270 


5706 


1.4 GHz continuum (mJy) 


10 


852 


<842 


4.8 GHz continuum (mJy) 


11-27 


1000 


<990 


8.3 GHz continuum (mJy) 


25 


800 


<775 


99 GHz continuum (mJy) 


23-42 


75 


<52 


Derived parameter 


Component I 




Component II 


Density (cm -3 ) 


~10 4 




300-1000 


Size (pc) 


<10 




30-60 


N Lyc (xltfV 1 ) 


4 




3-6 



Col (1): Parameter; col (2): predicted value of the parameter, summed over the eight compact H92a 
sources (Component I, see Sect. f£"2*l ; col (3): observed value of parameter, from low resolution data; col 
(4): the difference between col. 3 and col. 2, which is modelled as Components II and III in Sect. 1^31 



uum spectral index flatter, thus explaining the presence of 
thermal sources not detected in RRLs in our observations. 



6.1. The IR peak 

The IR peak dominates the NIR and MIR continuum, Br7 
and [Nell] lines and accounts for 10 % of the total bolo- 
metric IRAS flux from NGC 253, and is hence a prominent 
thermal source. However, this source is not conspicuous in 
the radio continuum image which is understandable since 
the total radio continuum is predo minantly non-thermal. 
Br7 data from lForbes et ahl l|l993j) implies a free-free flux 
density in the radio of only ~1 mJy in the for the IR peak 
for opti cally thin gas (and 2.5 mJy for A y=10); observa- 
tions of iKalas fc Wvnn- Williams] l|l994l) are consistent if 
a large MIR excess is present. However, the IR peak does 
not have a strong RRL counterpart, which is surprising. 
Much higher extinction towards the central 15 pc region 
than towards the IR peak could be invoked. Nevertheless, 
since the 12.8 /im [Nell] emission is much stronger towards 
the IR peak than towards the centre, the required extinc- 
tion towards the latter would need to be unrealistically 
high, unless most of the Ne in the centre is doubly ionised 
due to the AGN. However, low density models can be con- 
structed which can explain the weak RRL emission along 
with a high star formation rate. Hence we suggest that 
the thermal gas in the IR peak is probably at densities 
lower than ~few 100 cm" 3 and is therefore not promi- 
nent in high frequency line observations. We predict that 
a high resolution observation of the H166a line using the 
VLA in the A configuration would show the IR peak to 
be a prominent line source since this line is more sensitive 
to low density gas and the VLA in the A configuration 
will be able to resolve the IR peak from the surrounding 
emission. 



6.2. lonisation sources of cm-RRL emitting gas 

It was argued bv lMohan et alJ l)2002j) that the RRL emit- 
ting gas in the radio nucleus is ionised by a weak cen- 
tra l AGN, for which a dditional evidence was presented 
by I Weaver et alJ l|2002h . The nature of ionisation of the 
remainder of H92a emission is discussed in this section. 

If the individual H92a RRL sources in the 15 pc region 
are ionised by stars, the derived values of the ionisation 
rates and sizes imply the presence of a Super Star Cluster 
(SSC) within each (based on the typical sizes and lumi- 
nosities of SSCs derived bv lMeurer et alJll995|) . Similarly, 
massive clusters or SSCs would be required for the four 
sources outside the central region. However, the radio re- 
combination line widths (presumably caused by the ex- 
pansion of the HII regions) are supersonic, ^200 km s _1 
for the central sources and 70-100 km s _1 for the lat- 
ter ones. Hence the dynamical age for the ionised gas is 
~10 4 years for all of the sources in the central 15 pc re- 
gion and between (4-10) xlO 4 years for the extranuclear 
sources. These are too young to be realistic since it is im- 
probable that star clusters much younger than 10 5 years 
would have formed enough OB stars to provide an ion- 
isation in excess of 10 51 photons s — 1 . The thermal gas 
in the four nuclear sources is not gravitationally bound to 
the clusters within, though gravity could play a role in the 
dynamics of the gas in the extranuclear sources. This "life- 
time pro blem" is similar to th e situation in NGC 5253 and 
He 2-10 l|Mohan et all 1200 1|) . The actual age of the ob- 
jects could be increased by invoking internal mass loading 
whereby there is injection of mass into the ionised com- 
ponent from evaporation of neutral or molecular clumps. 
This can increase the densit y of gas within and maintain 
a constant Stromgren rad ius l|Dvson. William fc Redman! 
19951 iLizano et allll996|) . There could also be a substan- 
tial contribution from an unordered or turbulent compo- 
nent to the observed line width, in which case the dynam- 
ical time of the ionized gas will be longer compared to an 
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ordered outflow. However, since the central 15 pc region is 
also defi cient in PAH and MIR emission compared to the 
IR peak ijKeto et al.1 fl999") . the possibility of ionisation of 
these HII regions solely by SSCs is not entirely satisfac- 
tory. A possible explanation for the deficiency of nuclear 
PAH emission is dissociation by th e hard radiation from 
th e central AGN jKeto et all Il99flf) . 

I Weaver et al.l l|2002l) has shown that the radio nucleus 
is powered by a weakly accreting supermassive black hole 
or an intermediate mass black hole. Assuming the former, 
it is possible that the kinematics of the gas in the cen- 
tral 15 pes is influenced by the gravitational potential of 
the central mass. Hence a fraction of the large observed 
line widths could be explained by the AGN, though the 
remainder would still remain supersonic. The gas could 
also be partially ionised by the AGN, though star clus- 
ters would have to provide all of ionisation for the ex- 
tranuclear sources. There is evidence that the gas kine- 
matics in the central region of NG C 253 is complex (see 
iDas. Anantharamaiah fc YurJl200ll) . and involves multi- 
ple orbits at different orientations seen at high inclina- 
tion. Hence it is possible that the RRL sources in the 
15 pc region may not contain individual star clusters, but 
could be patterns or enhancements of emission measure 
along the line of sight. In such a case, the stellar den- 
sities needed for ionisation can be much lower, and the 
line widths of individual ionized gas components can also 
be smaller than the observed value, averting the "life- 
time problem" . However, the derived properties of the line 
emitting gas would remain unchanged. We predict that 
with increased spatial resolution, the line widths of the 
central RRL sources would decrease, similar to the [Nell] 
and RRL line widths in the Galactic Centre. In fact, due 
to the presence of a weak AGN, star clusters, high density 
ionised gas, and complex kinematics, the centre of NGC 
253 is qualitatively similar to the Galactic Centre. 

6.3. Multi-density model of the ionised gas 

Modelling the multi-frequency high resolution RRL data 
on NGC 253 led to a multi-densi ty model for the ionised 
gas, similar to that in Arp 220 ijAnantharamaiah et all 
2000). In contrast, using the low resolution RRL data 
at precisely the f requencies we have modelled here, 
IPuxlev et alJ l)l997[) could model all of the RRL emissions 
as arising from the same component of gas. The solution 
space of Components I and II in the n e -EM plane are plot- 
ted in Fig. El The figure for non-LTE conditions shows 
that the line spectral indices between 1.4 GHz and 99 GHz 
is not 2.0 for either component, though the spectral index 
of the total line emission is indeed very close to 2.0. Thus, 
high spatial resolution, coupled with multi- frequency data, 
is necessary to adequately model the compact ionised gas 
component. 

High resolution observations of molecular gas in 
NGC 253 within the central few hu ndred parse cs (HCN: 
IPaglione. Tosaki. Jacksor]ll995t CS: |Peng et aflll996|) in- 



dicate that the gas density exceeds 10 4 cm~ 3 . Hence 
Component I of ionised gas could well be produced by 
ionising this surrounding molecular material. Since the 
ionised gas becomes less dense and also larger in size with 
age, due to expansion of the HII regions, it is consistent 
to expect the more compact and denser Component I to 
have similar density as the molecular gas. The spectral in- 
dex of the global radio continuum emission in the central 

25 pc of NGC 253 is not steep ( 0.5; UA97, this work). 

From the spectral index of the continuum emission and 
also based on the [Nell] emission, we can conclude that 
free-free emission is not, however, the dominant portion 
of the total continuum emission. Therefore, the presence 
of foreground ionised gas with large covering factor (which 
absorbs the background non-thermal emission) is implied. 
Component II does satisfy these requirements, since it has 
a large transverse size and has sufficient optical depth at 
8 GHz. However, the exact location of this component of 
ionised gas needs to be investigated with future VLA A 
configuration observations of the 1.4 GHz RRL. Though 
the properties of Component III have not been derived 
here, this component can be modelled in detail using the 
observations of the 45 GHz RRL and also mm- wave RRLs, 
which have been carried out using the VLA and IRAM re- 
spectively (Zhao, and Viallefond, private communication). 

In this work, a three-component model of the ionised 
gas has been derived using observations of RRLs at three 
different frequencies. This raises the question of unique- 
ness of the decomposition into these specific density com- 
ponents. The ionised gas in NGC 253 is not expected to 
exist at three discrete densities; we interpret these results 
as an approximation to the intrinsic density distribution 
of the gas. This intrinsic distribution, presumably continu- 
ous, is a convolution of the three effects, namely, the initial 
density structure of the molecular ISM, the detailed time 
dependence of the star formation rate, and the dynamics 
of the HII regions. By observing widely spaced RRL, each 
of which is sensitive to a different range of densities, the 
properties of the ionised gas can be approximated in terms 
of a finite number of discrete density ranges. Further, the 
sizes and masses of the gas in each of these density ranges 
can also be predicted. Hence, we divide the ionised gas into 
three density bins, the exact ranges of which are decided 
partly by the frequencies of the observed RRLs. We can 
then accurately determine the length scales and masses 
and hence the required ionisation rates for the gas in each 
of these bins by modelling the line formation. Future ob- 
servations of a larger number of RRLs from decimetre to 
millimetre may be used to constrain a more general den- 
sity distribution function directly. 

7. Summary and Conclusions 

The 1.4 GHz line emission is detected with a resolution of 
-50 pc. The 8.3 GHz H92a and 15 GHz H75a line emis- 
sions have also been imaged, at high resolutions (—4 pc 
and ~2 pc respectively). The total H92a emission, which 
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is about half the value measured at low resolution, consists 
of an extended region of emission surrounding the radio 
nucleus (the '15 pc region') and four compact sources in 
the outer disk of NGC 253. The line emission in the 15 pc 
region can be decomposed into four additional compact 
sources. The peak of the line emission is towards the ra- 
dio nucleus. The line widths of the outer sources is <100 
km s _1 whereas the line widths of the sources in the in- 
ner 15 pc are about 200 km s" 1 . H75a emission has been 
detected from only two of these sources. Not all the H92a 
line sources have identifiable continuum counterparts. 

In contrast to existing interpretations, we show that 
multi-wavelength RRL emission from NGC 253 does not 
arise in optically thin gas in LTE. The four cxtranuclear 
H92a sources as well as the four compact sources in the 
15 pc region are modelled individually and arise in gas 
at densities close to 10 4 cm~ 3 . The total ionisation and 
gas mass in this high density Component I are ^4xl0 52 
s" 1 and 5000-10000 M Q . The total size subtended by 
Component I is <10 pc. This component produces a neg- 
ligible fraction of the 1.4 GHz line emission. 

A second component of ionised gas, Component II, 
is required to explain all of the observed H166a line at 
1.4 GHz and the remainder of the cm- wave RRL emis- 
sions. This gas is a lower density of a few 100 cm -3 . The 
The typical transverse size is a few tens of pc. A third 
component of gas, whose density must be higher than 10 4 
cm~ 3 , is also needed to explain the observed H40a emis- 
sion at 99 GHz, but has not been modelled in this work. 
In conclusion, 

— The ionised gas component in NGC 253 is decomposed 
into three distinct components: a low density compo- 
nent (~500 cm -3 ) which contributes most of the ob- 
served thermal radio continuum emission and almost 
all of the 1.4 GHz RRL; a higher density component 
(~10 4 cm' 3 ) which is more localised and exhibits su- 
personic line widths, and a third component at an even 
much higher density, responsible for mm- wave RRLs. 

— The nature and location of the IR peak is uncertain 
and probably contains low density gas ionised by a star 
cluster which is much less obscured than the central 
15 pc region. 

— The ionisation of the line emitting high density ther- 
mal gas in the 15 pc region is probably determined 
by obscured star clusters and the dynamics could 
be partly determined by the AGN, like the Galactic 
Centre. Alternatively, some of the compact line sources 
could be line of sight enhancements in emission mea- 
sure due to complex dynamics in the region, rather 
than individual HII regions surrounding massive star 
clusters. In either case, we have discovered further 
examples of supernebulae with large line widths sus- 
tained by internal gas dynamical processes. 
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Appendix A: Radio continuum sources in the disk 
of NGC 253 

The 1.4 GHz H166a observations performed us- 
ing the VLA in the C con figuration, described in 
lAnantharamaiah fc Gossl l)l990|) . were reanalysed. The 1.4 
GHz RRL was detected, and was found to coincide with 
the position of the peak of the radio continuum (radio 
nucleus). This dataset is used to image the continuum 
emission from NGC 253 at a higher sensitivity and over 
a larger field of view than what is available in the litera- 
ture. Sources within 20' of the centre of NGC 253 which 
are withi n the spatial extent of the 330 M Hz continuum 
emission l|Carilli et all Il992t ICarilll Il996|) are identified 
with the galaxy. Th ese are named a s SI to S9 and are 
marked in Fig. |A~TL luivestadl (2000) listed 22 circumnu- 
clear sources within the inner 3' (<2 kpc) region. Using 
our C as well as BnA configuration images, we mark those 
of their sources which have uncertain parameters, as well 
as new sources, in Fig. lA.lb . The flux densities and po- 
sitions of these continuum sources were determined us- 
ing JMFIT in AIPS and the catalogue list is presented in 
Table IA~T1 
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Table A.l. Continuum sources in the 1.4 GHz VLA C and BnA configuration images 



Source 


Position (B1950.0) 




£51.4 GHz 


Source 


Comments 


name 


a 


(00 h ) 


<* ( 


-25°) 


(mJy) 


size 




SI 


44 m 


32 s . 65 


40' 


11" 


.2 


3.4 


U° 




S2 


44 m 


35 s . 95 


26' 


52" 


.6 


10.8 


4.0 x 1.6 




S3 b 


44 m 


44 s . 57 


34' 


06" 


.0 


19.0 


0.5" x 0.4" 


Bl 


S4 6 


45 m 


17 s . 56 


30' 


00" 


.0 


25.4 




B2; 3 point sources 


S5 b 


45 m 


27 s . 50 


28' 


53" 


.0 


1.8 


U 


real feature ? 


S6 


45 m 


31 s .57 


34' 


53" 


.7 


56 


2.4" x 0.5" 


B3 


S7 


45 m 


34 s . 03 


43' 


59" 


.7 


7.2 


3"xl"? 


possibly resolved 


S8 


45 m 


33 s . 79 


32' 


49" 


.1 


1.6 


U 




S9 


45 m 


37 s . 70 


27' 


34" 


.9 


2.8 


U 




15° 


45 m 


06 s . 12 


31' 


55" 


.3 


2.0 


2"x2" 


possibly resolved 


19 c 


45 m 


12 s . 07 


32' 


55" 


.7 


1.2 


4.5"x2.l" 




20 c 


45 m 


13 s . 13 


33' 


01" 


.1 


3.5 


4.1"x3.4" 




22 c 


45 m 


15 s . 75 


32' 


01" 


.2 


1.5 


8.6" x 5.6" 


within diffuse em. 


23 d 


45 m 


12 s . 91 


31' 


23" 


.9 


0.6 




too weak 




45 m 


15 s .69 


31' 


51" 


.2 


0.9 




too weak 


25 d 


45 m 


16 s .80 


32' 


13" 


.9 


0.7 




too weak 



a U - unresolved, even in the BnA configuration images. 

^ These are within 2-3 synthesized beams of — ve extended emission in the C configuration image. 
c Sources in U00 lUlvcsiad 2000) with inaccurate flux densities, positions. 
Sources seen in BnA image but not in lUlvestadl 120001) . 
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Fig. A.l. (a) The 1.4 GHz continuum emission from 
NGC 253 taken using the VLA C configuration data, 
showing sources identified with the galaxy. The synthe- 
sized beam is 35"xl7" at a P. A. of 26°. The contour levels 
are in steps of (-20, -10, -7, 7, 10, 20, 40, 60, 80, 100, 200, 
400, 600, 800, 1000, 2000, 4000, 6000, 8000) times the rms, 
0.2 mJy/beam. (b) Higher resolution BnA con figuration 
image . The area shows the sources detected bv lUlvestadl 
l)2000() . The synthesized beam is 3.8"x2.6" at a P.A. of 
61°. The rms of the image is 0.1 mJy/beam and the con- 
tour levels are in steps of (5, 7, 10, 20, 40, 60, 80, 100, 
200, 400, 600, 800, 1000, 2000, 3000, 4000) times the rms. 
Negative contours are below — 5cr. 



